Population imbalance and pairing in the BCS-BEC crossover of three-component 

ultracold fermions 
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We investigate the phase diagram and the BCS-BEC crossover of a homogeneous three-component 
ultracold Fermi gas with a U(3) invariant attractive interaction. We show that the system at 
sufficiently low temperatures exhibits population imbalance, as well as fermionic pairing. We de- 
scribe the crossover in this system, connecting the weakly interacting BCS regime of the partially 
population-imbalanced fermion pairing state and the BEC limit with three weakly interacting species 
of molecules, including pairing fluctuations within a t-matrix calculation of the particle self-energies. 
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I. INTRODUCTION 

Multi-component ultracold atomic systems have re- 
cently been the focus of both experiment and the- 
ory, motivated in part by the prospect of simulat- 
ing a wider range of many-body models, such as lat- 
tice SU(N) models Jl|-|4j and quantum chromodynamics 
(QCD) analogs [|| ||, than is possible with single- or two- 
component systems. The possibility of creating analogs 
of color supcrfluid states and the formation of hadronic 
states in multicomponent systems [^-0 is especially in- 
teresting since the regime of cold dense QCD matter is 
not directly achievable in current nuclear experiments or 
in lattice QCD. 

When three species of fermions weakly attract each 
other, two species form Cooper pairs and the third re- 
mains a Fermi liquid § |l(J. Which two species pair de- 
pends on anisotropies in the interactions and mass differ- 
ences between different species. If there is no anisotropy, 
the Hamiltonian of the system possesses global U(3) sym- 
metry with respect to rotation in species space, and the 
pairing breaks this symmetry. An important feature 
of the three-component fermion system is spontaneous 
population imbalance, hrst noted in the continuum in 
Ref. [15 at T = 0. In addition, BCS superfluidity and 
population imbalance (magnetism), with two indepen- 
dent order parameters, can coexist, an intrinsic feature 
of a multicomponent Fermi systems, as shown by Cherng 
et al. in the weak-coupling BCS regime. 

We consider here U(3) invariant three-component ul- 
tracold Fermi gases in three-dimensional free space with 
varying interaction, with a focus on spontaneous popu- 
lation imbalance and superfluidity at finite temperature 
at general interaction strength; we study the phase dia- 
gram in general and the BCS-BEC crossover of the sys- 
tem, fixing only the total number of particles and allow- 
ing spontaneous population imbalance to occur. With 
a fixed total number of particles, population imbalance 
is accompanied by spatial inhomogeneities, such as, for 
example, domain formation. We first analyze the system 
at zero temperature in BCS mean field to show that the 
fermion pairing gap and population imbalance both de- 
velop with increasing bare attractive interaction between 



the fermions. Then we discuss nonzero temperature, 
starting from the BCS region where the scattering length 
is small and negative. We calculate the population im- 
balance as well as the BCS transition temperature there 
as a function of interaction strength and temperature, 
to lowest order in the interaction. The thermodynamic 
potential derived here agrees with previous calculations 
p"(j| n when the chemical potentials of the three species 
are equal. We also derive the Ginzburg-Landau free en- 
ergy as a function of the two order parameter-the pairing 
gap and the population imbalanc-and discuss a possi- 
ble analogy between dense QCD and three-component 
ultracold fermions. We then turn to the BEC limit of 
three-component ultracold fermions, where the scattering 
length is small and positive, a regime described by three 
different weakly interacting species of molecules made of 
different combinations of fermions. We show that Bose 
condensation of the molecules is accompanied by pop- 
ulation imbalance. Finally, we discuss the BCS-BEC 
crossover connecting BCS and BEC limits, following the 
procedure of Nozieres and Schmitt-Rink m7j to include 
pairing fluctuations (or non-condensed pairs), here in a 
summation of ladder diagrams for the self-energies; this 
calculation yields a transition temperature to the con- 
densate phase that reduces to the BCS and BEC limits. 

Degenerate three-component gases have been experi- 
mentally realized using the three lowest hyperfine states 
of 6 Li jl^, [l9|; at high magnetic fields, well beyond uni- 
tarity, the scattering lengths between the three hyperfine 
states are negative and sufficiently close that the sys- 
tem is approximately U(3) invariant. In addition, ul- 
tracold gases of alkaline-earth-metal atoms possess good 
SU(N) invariance (with N up to 10) and are good 

candidates to observe the physics discussed here. Ytter- 
bium has an SU(6) symmetry due to the nuclear spin; an 
SU(3) invariant mixture can be obtained by using only 
three spin components. In 6 Li as well as in m Yb and 
173 Yb, the temperatures currently achieved experimen- 
tally are around T > 0.3T F j§ [l9j. With a factor 
of ~ 3 decrease in temperature, phase separation due to 
the formation of population-imbalanced domains could 
be observed. 

Around the unitarity point, 1/a = 0, in a U(3) in- 
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variant system (where a is the s-wave scattering length) , 
three-body Efimov bound states can exist [po|-p7|. Efi- 
mov states have been experimentally observed in a trap 
through an increase of the particle loss rate, mediated by 
these states (l^, ^8|. In this paper, we analyze the 
system on time scales long enough to see the two-body 
interaction physics but short enough that Efimov states 
or three-body collisions can be neglected; such an inter- 
mediate thermalized regime can exist in a trap at suf- 
ficiently low densities, since the two-body collision rate 
is proportional to the particle density squared whereas 
the three-body collision rate is proportional to the den- 
sity cubed |5(J. As we show, the homogeneous state is 
unstable against the formation of inhomogeneous struc- 
tures with population imbalance; population imbalance 
suppresses the formation of Efimov states, tending to sta- 
bilize the inhomogeneous three-component system. 



II. THREE-COMPONENT U(3) INVARIANT 
FERMIONS 

We consider a three-component fermion system in free 
space with equal masses and the same scattering length 
between different species. We label the three species by 
"colors" in analogy with QCD, "red (r)," "green (g)," 
and "blue (b)." At low temperature, the interaction is 
dominated by s-wave scattering, and the Hamiltonian is 

k.a 
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where tp^ k is the creation operator of a particle with 
color a = r, g, b with momentum k; V is the volume, and 
we take h = 1 throughout. We assume an attractive bare 
contact interaction of strength U < 0. Although we take 
a common chemical potential /i for all three species, the 
numbers of each species in the state of lowest free energy 
can be different as a consequence of interactions, an effect 
that would be observable, in an experiment that starts 
with equal numbers, as an inhomogeneous state. The 
Hamiltonian is invariant under global U(3) rotations of 
the species. 

The attractive interaction leads to pairing of fermions 
at low temperature. The pairing order parameter is an- 
tisymmetric in color, and thus has the form 

A a (r) oc e ai a 7 (^>g(r)0 7 (r)}. (2) 

Since under a global U(3) rotation, 

ipa,k Uapipp^, (3) 

where U a p € U(3) (we use the convention that repeated 
indices are summed over), A a transforms as 

A Q (r) oc ^(^(r)^ 7 (r)) -> (det U)U* A^(r), (4) 



where ip a (r) is the Fourier transform of tp a ,k- 

To prove Eq. (|]), we consider the operator A Q = 
t a 1 '4>j3'4> 11 whose expectation value is proportional to 
A Q . The combination ip a A a transforms as 

ip T A = ip a A a = e a /3~ ( ?p a 4>i3?P~f 

= det U e^ip^ip v ip^ = det Ue a p^ip a ippip^ 

= det Uij?A. (5) 

On the other hand, ip T ^P T U T . Therefore, A — s- 
deti7(f7 T )- 1 = detUU*. 

As a consequence of the transformation (|J), we can 

- when the order parameter is independent of position 

- always choose appropriate axes of colors to transform 
the pairing order parameter into the form A = (0, 0, A), 
that is, by taking appropriate linear combinations of 
the species, we find that only two colors are paired and 
one is left unpaired. By applying a Bogoliubov-Valatin 
transformation, we can see that there are two gapped 
fcrmionic excitations corresponding to the quasiparticlcs 
of the paired fermions, and one ungapped excitation due 
to the unpaired fermions. In the following, we assume, 
without loss of generality, that the red and green particles 
are paired and the blue are not paired. 



III. BCS MEAN FIELD AT T = 

In this section, we consider the ground state of the 
system within mean-field BCS theory. We describe the 
pairing between r and g particles and unpaired b particles 
with the BCS-like ansatz, 

l*) = n( u k + «k< k ^,-k) II < k |vac), (6) 

k lkl<4 

where |wk| 2 + |fk| 2 = 1 ancl k h F is the b Fermi momentum. 
The parameters itk and i>k are determined by minimiz- 
ing (if?\7i — fjj\f\i£) at fixed fi. Following the standard 
procedure, we obtain 



6. 



(7) 



where ^k = k 2 /2m— \i and the gap A = — (U/V) itk^k 
is determined by 



A 



(8) 



We use the relation of the bare coupling U and the scat- 
tering length a |1 §01, 
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U ~ 4^ ~ V ^ fc 5 ' 



(9) 
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to rewrite the gap equation for A ^ in terms of a as 




1 



1 



(10) 



The chemical potential is determined by hxing the total 
number of particles N: 



N 



b \3 



k 



1 - 



V 



6tt 2 



(11) 



The same gap and number equations were derived in 
Ref. pl| using path-integral techniques. We solve the 
gap equation ([!(]) and the number equation ( |IT| ) simul- 
taneously to calculate the pairing gap and the number 
imbalance in terms of the scattering length. 

In Fig. Jy, we plot the pairing gap A, measured in units 
of tF — kp /2m, and the number of r particles N r divided 
by the total number of particles N, against — l/kpa, 
where kp = (67r 2 A r /3y) 1 / 3 . The right side of the figure 




FIG. 1: (Color online) The number of red particles divided by 
the number of total particles N r /N and the pairing gap A, in 
units of €f, at zero temperature, vs. —1/kpa. The solid line 
shows N r /N (left vertical axis) and the dotted line A (right 
vertical axis). 

corresponds to the weak-coupling regime (BCS region); 
the bare coupling becomes stronger toward the left side 
(BEC region) of the figure. As we see, |A| and the frac- 
tion of red particles, N r /N, increase with stronger inter- 
action. The N r /N axis ranges from 1/3 to 1/2; when 
N r /N — 1/3, all three species are equally populated, but 
for N r /N — 1/2, only r and g particles are present. In 
general, N r /N is greater than 1/3 in the interacting sys- 
tem, and it approaches 1/2 as the interaction becomes 
stronger. Thus the ground state of the interacting sys- 
tem always exhibits population imbalance, or magneti- 
zation (in analogy with a spin system). The magnetiza- 
tion arises physically through the gain of pairing energy 



when there are more particles in r and g states, and as 
remarked earlier, it would reveal itself in experiment as 
an inhomogeneous distribution of particle numbers. 

With this basic picture in mind, we turn now to 
nonzero temperature. 



IV. BCS REGION 

In the BCS region, where the scattering length a is 
negative and small, perturbation theory in terms of the 
scattering length well describes the system. We first de- 
rive the phase diagram in this region, and then we derive 
the corresponding Ginzburg-Landau free energy. 



A. Mean-field phase diagram 

The mean-field Hamiltonian T-Lm is 
H M ~ £«jV 

= E (& + if - N S) < k vvk 

- A* E ^nk^.-k - A E <-k<k 



Y.\A\ 2 - ^ (N r N g + N g N b + N b N r ) , (12) 



where 



U 



A = -— E(VV,k^g,-k)- 



(13) 



As was done earlier, we assume equal numbers of red 
and green particles, N r = N g . Also, we now include the 
Hartree energy, Uh = 47ra/m. Defining 



Uh 

6,k = 6c + ^f 2iv r , 



(14) 
(15) 



we rewrite the mean-field Hamiltonian as 
Um - fiAf 

= E (^r.k^r^k + ^feV'ff.k) + E 6>,kVikV">,k 

k k 

- A* E Vv,k^-k - A E <-k€,k 



V 

u 



&\ 2 -^{N r +2N r N b ), 



(16) 



which is essentially the BCS mean-field Hamiltonian for 
paired red and green particles plus normal blue particles. 
Diagonalizing by a Bogoliubov-Valatin transformation, 
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we find the thermodynamic potential 



,-Kb 



- £ (e k - fr.k) - y t \ A| 2 - ^ + 2iV r iV 6 ) , (17) 



where e fc = x /£ 2 k + |A| 2 . The condition 0O/d|A| 



gives the gap equation 

1_ v l-2/(e fc ) 



2£fc 



1 



(18) 



shaded region, A / and N r /N > 1/3, a small num- 
ber imbalance. We show in the next subsection using 
the Ginzburg-Landau free energy that A ^ implies 
N r /N > 1/3 and vice versa. To extend the theory to 
the unitarity and BEC regimes, we take pair fluctuations 
into account Jsif , in Sec. VI. 

In the next subsection, we derive the Ginzburg-Landau 
free energy of the system in the BCS regime, and derive 
the relations between the pairing gap and the number 
imbalance. 



B. Ginzburg-Landau free energy 



where f(x) — l/(e l3x + 1) is the Fermi distribution func- 
tion. Again, /j, is determined by the number equations 



tanh 



with 



N = 2N r + N b . 



pe k /2 \ 



£k 



(19) 
(20) 

(21) 



Numerically solving the gap equation (^) with the num- 
ber equation (pl|), we obtain the gap and number imbal- 
ance at given temperature and scattering length, shown 
in Fig. |[ The figure plots N r /N as a function of 




FIG. 2: (Color online) Phase diagram of the BCS region: 
N r /N vs. —1/kpa and temperature, in units of ef- The z- 
axis ranges from 1/3 to 0.35. The intersection of the surface 
and the bottom plane toward higher T is the transition line 
between the ordered and normal phases. 

— 1/kpa and T. The normal phase is the unshaded re- 
gion at higher T; here A = and N r /N = 1/3. In the 



The interplay between pairing and number imbalance 
is most easily seen from the Ginzburg-Landau free en- 
ergy, the expansion of the free energy in terms of the 
corresponding order parameters around the transition 
temperature. We define the order parameter for num- 
ber imbalance, 0, by 



N r 



N 



Fixing the total number of particles N 
have equivalently 



(22) 



2N r + Nb, we 



1 (Nb 
2\V 



N 
W 



(23) 



To derive the Ginzburg-Landau free energy it is conve- 
nient (in the derivation only) to let the chemical potential 
lib for b be different from the chemical potential /x r for 
r and g. The thermodynamic potential Cl(T, fj, r , fib) can 
be derived as in the previous subsection. The Helmholtz 
free energy is then 



F(A, 



Q + 2fj, r N r + n b N b , 



(24) 



in terms of which the Ginzburg-Landau (GL) free-energy 
density can be obtained by expanding 

J- GL (A,0) = i(F(A,0)-^(O,O)). (25) 



We define 



Ir.k = 77- — Mr + -77- (N r + N b ) (26) 
2m V 



~ k 2 U H 



(27) 



and ik = y £, 2 k + A 2 , and the chemical potential of the 
normal phase /iq implicitly through 



TV 

y 



1 



(28) 
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where f£ = k 2 /2m - ^ + 2U H N/3V. In terms of these 
quantities, the Ginzburg-Landau free-energy density is 



•Fgl(A,. 



-/3?i,,k 



(3V 



-«5 



v 

N 



+ 3U B <r + 2(/x r - + ^7 (2/*r + M6 - 3^o) 5 (29) 

in the expansion in <fi and A, we keep in mind that fi r and 
Hb are implicit functions of A and <p through the number 
equations 



V ^2 

k 



tanh /?£fe/2 



N 
W 



and 



V ^ e PSb.k + i 



AT 
37' 



^ (30) 



(31) 



The Ginzburg-Landau free energy up to fourth order in 
the order parameters is 



•Fgl(A, 



iA 2 



Cl 

C3^ 3 



A 4 + 3 — - U H 
■C40 4 -2— A 2 <j> + c 5 A 2 <p 2 , (32) 

Cl 



where c\ ~ C5 and 6 are all positive, but the sign of a 
depends on temperature. The detailed coefficients are 
given in Appendix A. 

The physically realized values of the order parameters 
minimize the Ginzburg-Landau energy; to leading order 
in the order parameters, we then have 



O0 \Cl 



GL 



dA 



2A 



a + 2 b 



-2^A 2 

Cl 

(c 2 ) 2 

Cl 



0, 



Co 

!— < 

Cl 



The first condition implies 



C-2 



3(1 - Cl U H ) 



A 2 



(33) 
0. (34) 

(35) 



indicating that if the pairing gap is nonzero, the num- 
ber imbalance is nonzero, and vice versa. The second 
condition, combined with Eq. (p5|) , implies 



A 



o + 2 6 + 2 



(c 2 ) 



,w 

ci 3ci(1-ciJ7h) 



A^ 



0. (36) 



In addition to the solution A = 0, when a < this 
equation has a second solution, with lower free energy, 



A 



2(6 + c 2 2 /ci - c 2 2 /(3ci(l - Cl U H ))) ' 



(37) 



The transition to fermion pairing is at the temperature 
at which a = 0. 

The Ginzburg-Landau free energy of three-component 
ultracold fermions has certain similarities to the 
Ginzburg-Landau free energy of dense QCD derived in 
Refs. p^-|35{, which makes multicomponent ultracold 
atoms a promising analog of dense QCD. The Ginzburg- 
Landau free energy of dense QCD has the form 



Cl Q cn(d,a) = ^\d\ 2 + ^\d\ 4 



C '„3 , b ' 4 



7 '|rf| 2 a + A'|d|V 



(38) 



where d is the quark-quark pairing order parameter and 
a is the chiral symmetry breaking order parameter. We 
attach primes to the coefficients to avoid possible confu- 
sion with similarly labeled quantities used earlier. The 
signs of a' and a' depend on the temperature and the 
strength of the couplings. As argued in Refs. p3|-|35[, /?', 
c', 7', and A' are positive. 

With the correspondence between the present system 
and the dense QCD system, A O d and <f> ■<-» a, we see 
that the two Ginzburg-Landau free energies have a sim- 
ilar structure. Although the original QCD Lagrangian 
has a local SU(3) gauge symmetry, the Ginzburg-Landau 
free energy (|38|), which does not take the gluonic de- 
grees of freedom explicitly into account, possesses only 
global SU(3) symmetry. To this extent, one can con- 
struct an analogy with ultracold atomic fermions. Simi- 
larly, Nambu-Jona-Lasinio models of QCD [[36| |37j also 
have only global SU(3) symmetry. Differences between 
the QCD free energy and that of ultracold fermions are 
that the sign of a' becomes negative at low temperature 
whereas the coefficient of <j) 2 is always positive, and in ad- 
dition the coefficients of a 3 and <fi 3 are opposite in sign. 
These differences are due to the fact that the dense QCD 
system can undergo chiral symmetry breaking without 
quark-quark pairing, but the three-component ultracold 
fermion system, beginning with equal populations, can- 
not spontaneously develop local number imbalance with- 
out fermion pairing; with the symmetric interaction we 
are assuming, number imbalance arises from the gain of 
pairing energy with an increasing number of paired par- 
ticles. It would be interesting to see how the analogy can 
be sharpened in multi-component atomic systems where 
spontaneous number imbalance and fermion pairing oc- 
cur independently, for example, with increased numbers 
of species or with deviations from fully symmetric inter- 
actions. 



V. BEC LIMIT 

We turn now to the BEC limit, where the scatter- 
ing length between fermions is small and positive. We 
can regard the system here as a collection of three types 
of weakly interacting bound Bose molecules, each made 
of two fermions, which can be red-green, green-blue, or 
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blue-red. The molecules Bose-condense at sufficiently 
low temperature. The condensate of molecules can be 
reduced to a condensate of one type of molecule by ap- 
propriately choosing the color axes, as with pairing in the 
BCS regime. The condensate in the BEC limit is com- 
posed of the same two colors that are paired in the BCS 
limit. 

At high temperature, the system is not condensed, but 
is simply a gas of thermally excited molecules. Unlike 
in the condensate, one cannot exclude the existence of 
three types of thermally excited molecules. Whether the 
high-temperature system develops a number imbalance 
depends upon the intermolecular interactions. For the 
same type of molecules, the effective scattering length is 
0.6a [|8|, where a is the scattering length of the con- 
stituent fermions. Between different molecules, as we 
show later, the effective scattering length is still 0.6a. 
Thus, above the condensation temperature, the system is 
described by three kinds of thermally excited molecules 
with the same interaction between all molecules. As we 
show in Appendix B, the uncondensed Bose system does 
not develop a spontaneous number imbalance as long 
as the interaction between the same types of bosons is 
greater than half of the interaction between the different 
bosons. Thus the present system does not exhibit number 
imbalance above the condensate transition temperature. 

We have, therefore, the following picture of the BEC 
limit. At high temperature the system is a homoge- 
neous mixture of three types of molecules. The Bose- 
Einstcin condensation temperature is that of nonintcr- 
acting bosons of mass 2m and density N/6V, 

Tbe ^mc^(^) 2/3w0 - 137T - (39) 

Below TbeCj the system is a mixture of the condensate 
of one type of molecule and a cloud of thermal molecules 
of three types, which vanishes at T = 0. 

We now show that the scattering length between dif- 
ferent molecules is the same as that, 0.6a, between like 
molecules. The derivation of Ref. |gg of the scatter- 
ing length between similar molecules depended on the 
symmetry of the four-particle scattering wave function. 
Since, as we show, the wave function for scattering of dif- 
ferent molecules has the same symmetry, the arguments 
of Ref. HU lead to the same scattering length. We write 
the four-particle scattering wave function between sim- 
ilar molecules, for example, red-green on red-green, as 
^(ri, r 2 ; r 3 , r 4 ), where ri denotes the position of the 
red fermion of the first molecule, r 2 is the position of the 
green fermion of the first molecule, r 3 is the red fermion 
of the second molecule, and r 4 is the green fermion of the 
second molecule. The symmetries due to Fermi statistics 
are 

* s (ri,r 2 ;r 3 ,r 4 ) = -* s (r 3 , r 2 ; ri, r 4 ) 

= -* s (r 1 ,r 4 ;r 3 ,r 2 ). (40) 

On the other hand, scattering between different 
molecules, for example, red-green and red-blue shown in 




FIG. 3: (Color online) Two different molecules colliding. 

Fig. |^, described by the four-particle scattering wavefunc- 
tion 4' ( i(ri, r 2 ; r 3 , r 4 ) (where r 4 now denotes the position 
of the blue fermion) , has only a single symmetry due to 
Fermi statistics, 

* d (ri,r 2 ;r 3 ,r 4 ) = -* d (r 3 , r 2 ; ri, r 4 ). (41) 

However, for s-wave scattering, the wave function is sym- 
metric with respect to the interchange of molecules, so 
that 

* d (ri,r 2 ;r 3 ,r 4 ) = * d (r 3 , r 4 ; r 4 , r 2 ). (42) 

Conditions Q) and (Q imply that 

*d(ri,r 2 ;r 3 ,r 4 ) = -* d (ri, r 4 ; r 3 , r 2 ), (43) 

which is exactly the same symmetry that was present due 
to the exchange of green fermions in \I/ S . 

The Schrodinger equation in the two cases has one ap- 
parent difference, that is, the delta-function interaction 
between the green and blue fermions. However, the an- 
tisymmetry ( |43|) for exchange of green and blue fermions 
implies that the product of the green-blue potential and 
the wave function in the Schrodinger equation vanishes, 
so that the Schrodinger equation is the same as for identi- 
cal molecules, and the scattering length is also the same. 
This argument depends crucially on the two molecules 
having one color (here red) in common. 

VI. CROSSOVER THEORY 

The crossover, in a two-component system, from BCS 
pairing in the weak-coupling region to a BEC of weakly 
interacting molecules in the strong-coupling region is con- 
tinuous, as seen in experiment |59|-fn| and understood 
theoretically || g§-§|. A common feature of the- 
ories of the BCS-BEC crossover at nonzero temperature 
is the incorporation of pairing fluctuations, which allow 
thermally excited Cooper pairs to exist above the con- 
densate transition temperature. We now apply this idea 
to develop a theory of the crossover, at nonzero tem- 
perature, in the three-component system to connect the 
BCS and BEC regimes discussed earlier, and see that 



7 



the crossover is also continuous. fl52| We incorporate pair- 
ing fluctuations through a self-consistent summation of 
ladder diagrams, and then numerically solve for the tran- 
sition temperature between the condensate and noncon- 
densate phases. 



A. Self-consistent summation of ladder diagrams 

We construct the crossover theory in terms of the finite 
temperature normal and anomalous Green's functions: 

G a (r -r',t- t') = -i (T (Mr, Wl(r', t')) ) 
J"(r - r',t - f) = -i (T (Mr, M g (v' , 0)> , ( 44 ) 

where T denotes time ordering. We assume still that 
pairing takes place between r and g particles. The pair- 
ing gap is given in terms of the Fourier transform of 
F(v - r',t- if) by 



equations, we obtain 



A = -U 



f d 3 k „ , U f d 3 k v-^ 



(45) 



where k denotes (k, LUk)', the summation is over the 
fermionic Matsubara frequencies uj^ = mVk/P with odd 
integer v^. The Schwinger-Dyson equations for the 
Green's functions, illustrated in Fig. |[ are 



Qr(k) = (go(k)- 1 ~ Err(A)- 



Yi r g(k)Yig r (k) 

g (-k)- x - E OT (-fc) 
^ fc(A:) = e (fc)-i-s bb (fc)' 

J- (k) = — Eg r (fe) • {E r g(fe)5]g r (fc) 

+(0o(fc)- 1 - E„(fc))(a (-fc) _1 - E 9S (-fc))} _1 . 

(47) 

The main contribution to the off-diagonal self-energies is 
the gap: 

J J f rPk' t , 

E ^ (fc) = J J ftrf^^ = ^) = - A ' ( 48 ) 

where we assume without loss of generality that A is real. 
Then the r-particle self-energy, for example, is given by 

/d 3 q 1 v 
J7j^3jj^2l T rg(k,k;q)g g (-k + q) 

+T rb (k,k;q)g b (-k + q)}, (49) 

where T a p(k,k';q), is the two-particle t matrix for in- 
coming particles of color a with momenta k and /3 with 
— k + q, and outgoing with momenta k' and — k' + q, 
respectively; the uj q are bosonic Matsubara frequencies. 
The corresponding diagram is Fig. |[ Including the t 




FIG. 4: The Schwinger-Dyson equations for the normal and 
anomalous Green's functions. 



g r (k) = g (k) + g (k) (x rr (k)g r (k) 
P(k) = g (-k) (-x gr (k)g r (k) + e 99 



rg (k)P(k)), 
Ub\ - r.J-b\ f-v !b\a OA j- E ao (-k)J*(k)) , 

g b (k) = g (k) + g (k)E bb (k)g b (k), (46) 



where Qo(k) = — £k is the free-particle Green's func- 
tion, and are self-energies with an incoming a par- 
ticle and an outgoing j3 particle. Solving this system of 



-->©--> 



>l 



I I 

i r i 

I rg i 



_> — 



1 rb 



FIG. 5: Self-energy written in terms of t matrices. 

matrix in the self-energy takes pairing fluctuations into 
account, and as shown in Ref. |T7| ], encompasses thermal 
fluctuations of paired molecules in the BEC limit and the 
Hartree approximation in the BCS limit [jl^], thus con- 
necting both limits continuously. Note that there is no 
process of this form in which the top line is anomalous 
since such a process would involve scattering between two 
r particles, either initially or finally, which is forbidden 
by the Pauli principle; the internal lines can, however, be 
anomalous. 

On the other hand, in the self-energy of b particles, 
the top line can in principle be anomalous; however, this 
process would involve particle-hole scatterings either ini- 
tially or finally, which is negligible for short-range inter- 
actions J4!l; the self-energy involves only a sum of rb 



= -U-U 



and gb particle-particle scatterings. The Bethe-Salpeter 
equation for the rb t matrix becomes 

r rb (fc,fc';q) 

/ 7^^3 ^ XI + Q') r T'6(p, k'-q). 

(50) 

As one sees by iterating this equation, T rb (k, k'; q) is inde- 
pendent of k and k'\ we write T(k,k';q) — T(q). Solving 
Eq. (|50|), we obtain 

' (51) 

Tgb takes the same form mutatis mutandis. 

In T rg we must take the rg anomalous Green's func- 
tions into account, as illustrated in Fig. ^. Solving the 



g g r r g g 
> — i~> — ^-r< — >n — 



-5 — ■— s — — >-■ — > — 
r g g r r 



FIG. 6: An anomalous contribution to the rg £ matrix. 

Bethe-Salpeter equation in Nambu matrix notation, we 
have 



r , s Xu(-g) 

r9{q) Xii(9)Xii(-g)-Xi2(9) 2 ' 



(52) 



where 



Xi2(g) 



*7 J (2ir) 3 /3 
d 3 p 1 



(54) 



To determine the gap and the number imbalance as 
a function of temperature and scattering length involves 
self-consistently solving the gap equation ( [47|) , which can 
be rewritten as 

i _ r d 3 k 1 



(Go(k)- 1 - S rr (fc))(g (-fc)- 1 - £«(-*)) + A 2 ' 

(55) 

together with the number equations 

X> ialfc7? Sr(fe), (56) 



N r 


= lim / 


f d 3 fc 


1 


V 


(2tt) 3 


1 


N b 


= lim 1 


f d 3 k 


1 


V 


(2tt) 3 





However in this paper we focus only on calculating the 
transition temperature. 



B. Evaluation of T c 

We now use the formalism of the previous subsection 
to evaluate the transition temperature, where the pairing 
gap A becomes zero. The gap equation at T c is equivalent 
to the condition that T rg (q) diverges at q = 0. Therefore, 
at T c , we can make the approximations, 



(k) 



d 3 q 1 . 

TTT-^T- E ( V rg{q)Qg{q - k) + T rb {q)g b (q - k)) 

£ E ( T r g (q)G g (-k) + T rb (q)g b (-k)) , 

(58) 



and 



^bb(k) 



d 3 q 2 



(2^) 3 /?, 
d 3 q 2 



^r br (g)^ r .(g - fc) 

{27T) s pX™^- (59) 

For T > T c , the t- matrices do not depend on the color 
indices. Then, using the final line of Eq. (|58| ) we see that 
the Green's function for r particles becomes 

where we introduce a "pseudogap" A pg at T c by writing 



r/ 3 



9 \ - 



A J (2tt) 3 



r(«), 



(61) 



with r = r,. s = T rb = T br . The final line of Eq. ( p0| ) is 
just a BCS Green's function with the gap replaced by the 



pseudogap. We write -Ek = yCk + A 2 g for convenience. 
Similarly Q b (k) has the same form at T = T c . 
The number equations then reduce to 



W 2 J (2tt) 3 V £k 2 
while the equation for the pseudogap is 
1 f d 3 k 1 



(62) 



t/ J (2tt) 3 /3 ( 



d 3 fc 



(27T) 



E^( fc )^(- fc ) 

£ 2 \ tanh(/3 c £;k/2) A 2 g /'(£ k ) 



El 



El 2 



(63) 
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as before, the bare coupling U is related to the scattering 
length a through Eq. (||). 

In the BCS limit, kpa — > 0~, A 2 g tends to zero, as 
we can see by considering the BCS gap equation at T c 
(not the mean-field BCS transition temperature, but the 
same T c that we are using here) with a gap A 



0.25 



1 

U 



d 3 k jtanh(/3 cV /£j 



AV2) 



(20 



2Ve 



(64) 



Expanding the right sides of (|3j) and (|64| ) in terms of A 2 S 
and A 2 , we see that the zeroth order terms are identical. 
Also, since the final line of Eq. (p3[) decreases monoton- 



ically with A 2 g , the limit A — > 0, as in weak-coupling 
BCS, implies A 2 pg -> 0. 

Determining T c requires estimating A 2 g , which we do 
by expanding T rg (q)~ 1 around q = 0, recalling that 
rV^O)" 1 = at T c : 



1 



-rv 3 (q,w 9 ) = — 



d 3 p 1 



U ' J (2tt)3 p c 
d 3 p 1 



(2tt)3 A 



1 



k=—p 



= Zuj q - 79 . 



(65) 



Explicit forms for if and 7 are given in Appendix C. The 
pseudogap then becomes 



A 2 - -2 



1 



1 

d 3 g 1 _ C(3/2) 

Z 7 (2tt) 3 e^7<? 2 /2 - l 4Z 



7T/?c7 



3/2 



(66) 



Solving the number equation (p2|), the gap equation 
(|63|), and the expression for the pseudogap (|6^) self- 
consistcntly, we obtain the transition temperature, plot- 
ted against —1/kpa in Fig. 0. The solid line in the 
figure is the transition temperature calculated with the 
ladder summation formalism described here, and the dot- 
ted line is the result from mean-field BCS theory. The 
ladder summation line approaches the mean-field line in 
the BCS limit. On the other hand, in the BEC limit, 
the ladder summation correctly yields T c 0.137 Tp. 
The crossover theory presented here connects both limits 
continuously. 

Throughout, we have kept a common chemical poten- 
tial for the different species, and found that below T c the 
number of b particles becomes smaller than the number 
of r or g particles. In ultracold atomic experiments, the 
number of the particles in each species is usually fixed 
at the start, and thus the simplest scenario that may 




FIG. 7: (Color online) The phase diagram of three-component 
ultracold Fermi gas. The temperature is in units of Tf- The 
solid line is the transition temperature calculated with pair- 
ing fluctuations incorporated through the summation of lad- 
der diagrams. The dotted line is the transition temperature 
calculated from mean-field BCS theory. The mean-field line 
corresponds roughly to the temperature at which fermions 
start to form (noncondensed) pairs. The line calculated from 
the ladder summation is where the Cooper pairs start to con- 
dense. Toward the left end of the figure, the transition tem- 
perature approaches the BEC limiting value T c ~ 0.f37TF. 



occur experimentally is that the number imbalance ap- 
pears through the formation of population-imbalanced 
domains. The formation of population-imbalanced do- 
mains leads to a gain of condensation energy of order 
E c V/2 for the fully imbalanced state, where E c is the 
condensate energy density in a balanced system; the fac- 
tor 1/2 = 3/2 — 1 is the increase in the relative number 
of Cooper pairs in the imbalanced state over that in the 
balanced state. On the other hand, the formation of a 
single domain wall costs a net surface energy E sur f of 
order E C V£, C /L, where £ c is the coherence length and L 
is the linear size of the system. The condition that the 
formation of the domain is beneficial for the system is 
E c V/2 > E surf , or roughly L/£ c > 1, which typically 
holds well. Domain formation is expected to decrease 
the free energy from that of the homogeneous state at 
low temperature. Other possible realizations of popu- 
lation imbalance include the formation of a "color den- 
sity" wave or the formation of an inhomogeneous (Fulde- 
Ferrell-Larkin-Ovchinnikov) superfluid; we leave analysis 
of these states as future study. Also to apply the present 
theory quantitatively under realistic experimental con- 
ditions it will be necessary to investigate the effects of 
Efimov states. 
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Appendix A: The coefficients of the 
Ginzburg-Landau free energy 

We outline here the derivation of the Ginzburg-Landau 
free energy ( |32| ) from the free energy J 7 gl(A, </>), Eq. (p9[). 
Since A always appears squared in the equations, odd 
powers of A do not occur in the free energy. To find 
the coefficients of A 2 and A 4 , we set = 0, and expand 
.Fgl(A,0) in powers of A 2 . Taking the derivative of the 
number equation for blue particles (^lj) with respect to 
A 2 , we see that fib (here allowed to differ from fi r ) does 
not depend on A 2 . Differentiating the number equation 
for red particles (|30|), we obtain 



dfi r 



OA 2 



£2 

ci ' 



where the subscript denotes the derivative at A 
0, and 



C2 



1 E 

v 4^ 



Aanh/3^/2 , 



V 4(£«) 2 



Note that both c\ and C2 are positive. Then 
clF GL (A,0) 

' k 

(c 2 ) 2 



<9A 2 

1 9 2 J- GL (A,0) 
2 



1 1 tanh 



U V 



a(A 2 ) 2 



= b 



Cl 



where 



Aanh/3^/2 , f(£°) 



4(£ k ) 2 



> 0. 



We similarly derive the coefficients of 



9J"gl(0» 


= 0, 


d(t> 





1 5 2 J- GL (O,0) 




2 dcj> 2 


V c i 


1 5 3 J- G l(O,0) 


«1 


6 90 3 


(ci) 3 ' 


1 a 4 J- GL (O,0) 


3 


24 <9</> 4 


4( Cl ) 4 



>l) 2 



Cl 



«2 



(Al) 

= = 

(A2) 
(A3) 

(A4) 
(A5) 

(A6) 
3 , and 

(A7) 
(A8) 
(A9) 
(A10) 



where 

k k 

Finally, the coefficients of 0A 2 and (f> 2 A 2 are 
5 2 -F G l(A,<? 



d^dA 2 
1 9 3 J- GL (A,, 



A=rf>=0 



2 {d4>) 2 dA 2 

- C2Kl _ _ 
" (ci) 3 ( Cl 



A=A=0 



!) 2 



Co 

!— 

'ci ' 



2^ 



(All) 



(A12) 



46 



(A13) 



Therefore, the Ginzburg-Landau free energy up to fourth 
order in the order parameters is 

•F GL (A,</>) 



1 1 ^ tanh 0£/2 \ A2 



U V 



6+M')A 4 



ci 



+ ■ '> I - - U H 

Cl 



+ 



4(d) 



Cl 



2g 

«i 
" (ci) 2 



K2 U4 _ 2 £2 + c5 ^2 A 2 



= aA 2 +(6+M) A 4 
+ c 3 3 + c 4 4 -2-A 2 ^ 

Cl 



ci 

2>(--U H 

\ci 

c 5 A 2 </> 2 . 



(A14) 



Note that the Ci and b are all positive. Also, since Uh is 
negative, the coefficient of 2 is positive. 



Appendix B: Population imbalance in a Bose 
mixture above the condensation temperature 

We derive the condition for the homogeneous state 
with population balance to be stable. Although the 
three-component ultracold Fermi gas can form three 
types of molecules, the basic physics of the instability 
toward inhomogeneous states can be captured by consid- 
ering a two-component Bose system. 

We derive the Ginzburg-Landau free energy of a system 
of two species of bosons, a and b, as a function of their 
population imbalance at fixed total number N = N a +Nb- 
With ak and 6k the annihilation operators of bosons a 
and 6 of momentum k, the Hamiltonian is 



H - fi a N a - fi b N b 

'*! 
2m 



= £ ( t - a k fl k + £ ( - Mfe ) bib* 



£ (4+q a L-q a k'ak + 6j i+q 6j c ,_ q 6k'6k 



2V 

Y £ a k+q 6 k'-q & k'ak, 



k : k',q 



(Bl) 



k,k',q 



11 



where Uq = Airao/m and U\ = Altai /m are the s-wave 
interaction strength between the same type and between 
different types of bosons, and ao and a% are the corre- 
sponding scattering lengths. 

We assume a sufficiently high temperature that neither 
system is condensed. In the Hartree-Fock approximation, 
we obtain 

H - Li a N a - fi b N b 



E 



;&.k 



2U Q 



V 

Nb 
V 



V 



y(N 2 a +NZ)-^N a N b , 



(B2) 



where £ a! k = k 2 /2m — fi a and = k 2 /2m — [i b . The 
number of particles N a and N b satisfies the self-consistent 
equations, 

N a = 9 (&,k + 2U n a + U x n b ) (B3) 

k 

N b = Y,9 (6,k + 2U n b + Um a ) , (B4) 

k 

where g(x) = l/(e@ x — 1) is the Bose distribution func- 
tion, and n a — N a /V and n b — N b /V. Then the ther- 
modynamic potential is 

y = -U (n 2 a + n\) - Uin a n b 

+ Jy X) 111 i 1 ~ eX P (-P (6»,k + 2U n a + UWb))} 

+ -^y J2 In {1 - exp (-/3 (6, k + 2U n b + Um a ))} , 
1 k 

(B5) 

and the Helmholtz free energy is 



F n 

— = — + fi a n a + Hbn b . 



(B6) 



The condition for the stability of the homogeneous 
state is found by expanding the Helmholtz free energy 
in terms of the deviation of the number of particles from 
the homogeneous state. We write the deviation of the 
numbers of particles from the balanced case as 



<P = n a -- = 



n b 

2 



then 



d F 

d 2 F 
dip 2 V 



tp=0 



2 [2U - U t 



G 



(B8) 
(B9) 



where 



^iE^^-w+^ + ^l (BI0 » 

k ^ 



and the homogeneous chemical potential /io is determined 
by 



n 1 v-^ / k 2 n n 

9 = yI>(^-^ + 2 ^9 +C/l 9 



k 

2m 



(Bll) 



The homogeneous state is stable if and only if 
d 2 (F/V)/dip 2 > 0. Since G < 0, we immediately con- 
clude that when 2Uq > U\, as in the present system, the 
homogeneous state is always stable at T > Tbec- [For 
2Uo < Ui, one finds G —> 0~ as T — > oo, and G — » — oo as 
T approaches Tbec from above, implying a phase tran- 
sition from the homogeneous to an inhomogeneous state 
at T > Tbec- The transition temperature increases with 
increasing U\ — 2Uq. As Ui — > 2Uq from above, the tran- 
sition temperature approaches Tbec from above.] 

Since the interaction is the same as that between 
identical and different molecules in the BEC limit of 
three-component ultracold fermions, the result derived 
here implies that the system is homogeneous above the 
condensation temperature. 



Appendix C: Expansion of IY 9 (q, uj q ) 1 

The expansion of r rg (q, ujg)^ 1 can be explicitly car- 
ried out using Eq. (|60|), with the result of Eq. (65), 
-r r9 (q, LOqy 1 ~ Zio q - 7<j 2 , where 



d 3 k f tanh(/3 c g k /2) A Ck rpi s 

Z=l j2^{ 2K + /(i?k) )^2 (CI) 



2EI 



and 



d 3 k 



(2tt) 3 2mEl 



?k P9 3m 8 



£ k T 2 (2£ 2 -A 



2 ) 
pgJ 



( B7 ) +^ (Al 9 - g) (El + £) } (tanh ^ + 2E k f(E k ] 

+ { + l£ - + ^) | Etf(E k ) 



(C2) 



12 



[1] C. Honerkamp and W. Hofstetter, Phys. Rev. Lett. 92, 

170403 (2004); Phys. Rev. B 70, 094521 (2004). 
[2] M. A. Cazalilla, A. F. Ho, and M. Ueda, New J. Phys. 

11, 103033 (2009). 
[3] A. V. Gorshkov, M. Hermele, V. Gurarie, C. Xu, P. Juli- 
enne, J. Ye, P. Zoller, E. Demler, M. Lukin, and A. M. 

Rey, Nat. Phys. 6, 289 (2010). 
[4] S. Taie, Y. Takasu, S. Sugawa, R. Yamazaki, T. Tsuji- 

moto, R. Murakami, and Y. Takahashi, Phys. Rev. Lett. 

105 190401 (2010). 
[5] A. Rapp, G. Zarand, C. Honerkamp, and W. Hofstetter, 

Phys. Rev. Lett. 98, 160405 (2007). 
[6] A. Rapp, W. Hofstetter, and G. Zarand, Phys. Rev. B 

77, 144520 (2008). 
[7] F. Wilczek, Nat. Phys. 3, 375 (2007). 
[8] K. Maeda, G. Baym, and T. Hatsuda, Phys. Rev. Lett. 

103, 085301 (2009). 
[9] A. G. K. Modawi and A. J. Leggett, J. Low Temp. Phys. 

109, 625 (1997). 
[10] P. F. Bedaque and J. P. D'Incao, Ann. Phys. (NY) 324, 

1736 (2009). 

[11] T. Paananen, J. -P. Martikainen, and P. Torma, Phys. 

Rev. A 73, 053606 (2006). 
[12] T. Paananen, P. Torma, and J.-P. Martikainen, Phys. 

Rev. A 75, 023622 (2007). 
[13] H. Zhai, Phys. Rev. A 75, 031603 (2007). 
[14] G. Catelani and E. A. Yuzbashyan, Phys. Rev. A 78, 

033615 (2008). 

[15] L. He, M. Jin, and P. Zhuang, Phys. Rev. A 74, 033604 

(2006) . 

[16] R. W. Cherng, G. Refael, and E. Demler, Phys. Rev. 

Lett. 99, 130406 (2007). 
[17] P. Nozieres and S. Schmitt-Rink, J. Low Temp. Phys. 59, 

195 (1984). 

[18] T. B. Ottenstein, T. Lompe, M. Kohnen, A. N. Wenz, 

and S. Jochim, Phys. Rev. Lett. 101, 203202 (2008). 
[19] J. H. Huckans, J. R. Williams, E. L. Hazlett, R. W. Stites, 

and K. M. O'Hara, Phys. Rev. Lett. 102, 165302 (2009). 
[20] V. Efimov, Phys. Lett. B 33, 563 (1970); Nucl. Phys. A 

210, 157 (1973). 
[21] E. Braaten and H. Hammer, Phys. Rep. 428, 259 (2006); 

Ann. Phys. (NY) 322, 120 (2007). 
[22] S. Floerchinger, R. Schmidt, S. Moroz, and C. Wetterich, 

Phys. Rev. A 79, 013603 (2009). 
[23] S. Moroz, S. Floerchinger, R. Schmidt, and C. Wetterich, 

Phys. Rev. A 79, 042705 (2009). 
[24] S. Floerchinger, R. Schmidt, and C. Wetterich, Phys. 

Rev. A 79, 053633 (2009). 
[25] T. Luu and A. Schwenk, Phys. Rev. Lett. 98, 103202 

(2007) . 

[26] P. Naidon and M. Ueda, Phys. Rev. Lett. 103, 073203 
(2009). 

[27] E. Braaten, H. W. Hammer, D. Kang, and L. Platter, 

Phys. Rev. Lett. 103, 073202 (2009). 
[28] T. Kraemer, M. Mark, P. Waldburger, J. G. Danzl, C. 

Chin, B. Engeser, A. D. Lange, K. Pilch, A. Jaakkola, 

H.-C. Nagerl, and R. Grimm, Nature (London) 440, 315 

(2006). 

[29] V. M. Galitskii, Zh. Eksp. Teor. Fiz. 34, 151 (1958) [Sov. 



Phys.-JETP 7, 104 (1958)]. 
[30] A.J. Leggett, in Modern Trends in the Theory of Con- 
densed Matter, edited by A. Pekalski and R. Przystawa, 

(Springer- Verlag, Berlin, 1980). 
[31] J.-P. Martikainen, J. J. Kinnunen, P. Torma, and C. J. 

Pethick, Phys. Rev. Lett. 103, 260403 (2009). 
[32] M. Yu. Kagan and A. V. Chubukov, Pisfma Zh. Eksp. 

Teor. Fiz. 47, 525 (1988) [JETP Lett. 47, 614 (1988)]. 
[33] T. Hatsuda, M. Tachibana, N. Yamamoto, and G. Baym, 

Phys. Rev. Lett. 97, 122001 (2006). 
[34] N. Yamamoto, M. Tachibana, T. Hatsuda, and G. Baym, 

Phys. Rev. D 76, 074001 (2007). 
[35] G. Baym, T. Hatsuda, M. Tachibana, and N. Yamamoto, 

J. Phys. G 35, 104021 (2008). 
[36] T. Hatsuda and T. Kunihiro, Phys. Rep. 247, 221 (1994). 
[37] H. Abuki, G. Baym, T. Hatsuda, and N. Yamamoto, 

Phys. Rev. D 81, 125010 (2010). 
[38] D. S. Petrov, C. Salomon, and G. V. Shlyapnikov, Phys. 

Rev. Lett. 93, 090404 (2004). 
[39] M. Greiner, C. A. Regal, and D. S. Jin, Nature (London) 

426, 537 (2003). 
[40] M. W. Zwierlein, C. A. Stan, C. H. Schunck, S. M. F. 

Raupach, A. J. Kerman, and W. Ketterle, Phys. Rev. 

Lett. 92, 120403 (2004). 
[41] M. Bartenstein, A. Altmeyer, S. Riedl, S. Jochim, C. 

Chin, J. H. Denschlag, and R. Grimm, Phys. Rev. Lett. 

92, 120401 (2004). 
[42] R. Haussmann, Z. Phys. B 91, 291 (1992). 
[43] C. A. R. Sa de Melo, M. Randeria, and J. R. Engelbrecht, 

Phys. Rev. Lett. 71, 3202 (1993). 
[44] P. Pieri, L. Pisani, and G. C. Strinati, Phys. Rev. B 70, 

094508 (2004). 

[45] H. Hu, P. D. Drummond, and X.-J. Liu, Nat. Phys. 3, 
469 (2007). 

[46] K. Levin, Q. Chen, C.-C. Chien, and Y. He, Ann. Phys. 
(NY) 325, 233 (2010). 

[47] R. Haussmann, Self- consistent Quantum-Field Theory 
and Bosonization for Strongly Correlated Electron Sys- 
tems (Springer, Berlin, 1999). 

[48] P. Pieri, L. Pisani, and G. C. Strinati, Phys. Rev. B 70, 
094508 (2004). 

[49] A. Fetter and J. D. Walecka, Quantum Theory of Many- 
Particle Systems (Dover, New York, 2003). 

[50] Huckans et al. |19| argue that strong interactions with 
a long lifetime (> 0.1s) can in fact be achieved in a low- 
density gas. 

[51] We have so far assumed that the blue particles do not 
pair. However, the blue particles feel an effective attrac- 
tive interaction with each other mediated by the exis- 
tence of the red and green particles which can lead 
to p-wave pairing state of the blueparticles. However, as 
shown by Kagan and Chubukov [B2|, the transition tem- 
perature to such p-wave pairing is too low (T c ~ 10 -7 2>) 
to be observed in experiment, and we ignore it here. 

[52] At sufficiently low temperature Efimov states can lower 
the energy around unitarity, producing a discontinuous 
transition from the BCS to the BEC regimes |22y. 



